Regge gravity from spinfoams 
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We consider spinfoam quantum gravity in the double scaling limit 7 — > 0, j — > 00 with -yj = const., 
where 7 is the Immirzi parameter, j is the spin and -yj gives the physical area in Planck units. We 
show how in this regime the partition function for a 2-complex takes the form of a path integral over 
continuous Regge metrics and enforces Einstein equations in the semiclassical regime. The Immirzi 
parameter must be considered as dynamical in the sense that it runs towards zero when the small 
wavelengths are integrated out. In addition to quantum corrections which vanish for h — » 0, we find 
new corrections due to the discreteness of geometric spectra which is controlled by 7. 
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INTRODUCTION 

Spinfoams [1-3] are a tentative covariant quantization 
of general relativity. They provide transition amplitudes 
between quantum states of 3-gcomctry, in the form of a 
Misner-Hawking sum over virtual geometries [4, 5]. In 
the so-called 'new models' [6-8], intermediate quantum 
states are the ones of canonical loop quantum gravity, 
the SU(2) spin-network states, a remarkable feature that 
promotes the spinfoam framework to a tentative path 
integral representation of loop quantum gravity. 

The physical picture emerging from the spinfoam grav- 
ity is the following: spacetime is a quantum foam of vir- 
tual geometries with a discrete and purely combinatorial, 
relational structure, where the Planck scale plays the role 
of a natural minimal length. Spinfoams are the result of 
the quantization of general relativity formulated as a con- 
strained BF theory. In the BF theory with gauge group 
G, the basic variables are a 2-form B and a connection 
1-form to, both valued in the Lie algebra of G. The equa- 
tions of motion of BF theory impose the flatness of the 
connection 



F(oj) =du> + u>Auj = 



(1) 



where F is the curvature. So this theory has no local de- 
grees of freedom. The simplicity constraints on the field 
B turn the topological BF theory into vacuum general 
relativity, described by the Hoist action 

S = J (e A e)* A F(u) + - J (e A e) A F(oj) (2) 

with gauge group 50(1,3) (or 5*0(4) in the Euclidean 
signature). A trace in (2) is understood and * is the 
Hodge dual in the internal Minkowski (or Euclidean) 
space. The first term is the Einstein-Hilbert-Palatini ac- 
tion for general relativity in terms of the cotetrad e and 
the connection u>, regarded as independent variables. The 
second term vanishes on-shell by using the Cartan struc- 
ture equation D^e = 0, obtained varying uj. Varying e, 



(3) 



we have the Einstein field equations 

ei.iK L e J A F(lu) kl = 0, 

equivalent to the Ricci-flatness i? M „ = of the metric. 

The real number 7 ^ is called Barbero-Immirzi, or 
Immirzi parameter, and controls the spacetime discrete- 
ness as it enters the discrete spectra of area and volume 
operators [9, 10]. Spinfoam models provide a Feynman 
path integral, or state sum, based on a discretization of 
(2) over a 2-complex (generalized triangulation) and the 
full degrees of freedom are recovered in the infinite re- 
finement limit or equivalently summing over 2-complexes 
[11]. The spinfoam theory is sufficiently simple [12] and 
possesses the correct symmetries [13]. Recently it has 
been successfully coupled to matter fields [14, 15]. Fur- 
thermore, the large distance analyses were able to extract 
the correct low-energy physics in some simple cases [16- 
19]. 

In spite of those encouraging results, we are still lacking 
some conceptual steps in the derivation of the covariant 
form of loop quantum gravity. However, the opposite 
point of view is to take the quantum theory that has 
been 'derived' in various ways, e.g. in the form [6-8], 
and check for its physical viability (see section IV of [20] 
on this viewpoint). A major open problem is to show 
that the spinfoam Feynman path integral contains Ein- 
stein equations in some semiclassical regime. What we 
would like to have is the spinfoam version of the following 
semiclassical expansion of the gravitational path integral 



,is BJS ( S °„) 



(4) 



in the classical limit h — > 0, where Seh is the Eintein- 
Hilbcrt action for general relativity and on the right hand 
side it is evaluated on the classical solution <7°„ of the 
equations of motion determined by the boundary con- 
ditions specified on the metric g^ v . Here we propose a 
solution to this problem. 

In a previous paper [21] the authors introduced a sim- 
ple triangulation of spacetime with three dual vertices 
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(three 4-simplices) and one internal face (dual to the 
hinge, where curvature is concentrated) and showed that 
the spinfoam boundary amplitude is peaked on geome- 
tries with nonzero curvature, but the result was limited 
to the geometry considered and relied strongly on a par- 
ticular, and somewhat artificial semiclassical limit. In 
this paper we improve and generalize the previous argu- 
ments to a general 2-complex. In the continuous spec- 
trum limit 7 — > and in the semiclassical limit h — > we 
find the analogous of the WKB expansion (4) for spin- 
foam quantum general relativity. The analysis is based 
on the path integral formulation of references [22, 23]. 

The paper is organized as follows. First, we review 
the definition of the EPRL/FK spinfoam model in the 
Bloch coherent state basis for angular momentum. In 
the subsequent two sections we perform the continuum 
and semiclassical limit. After discussing the extension to 
the Lorentzian signature, we generalize the analysis to 
the case where a boundary state is present, and motivate 
an asymptotic formula for the spinfoam transition ampli- 
tudes in the semiclassical regime. In the last section we 
give some outlooks for the future developments. 

Throughout the paper we work in natural units G = 
h = c = 1 , but will restore some of those constants when 
needed. In particular, restoring only h the Planck length 
is lp — \f%. 



THE SPINFOAM AMPLITUDE 

We consider the spinfoam amplitude [6, 8] for a 2- 
complex a without matter. We restrict for simplicity to 
the Euclidean signature and to Barbero-Immirzi param- 
eter < 7 < 1, where formulas get simpler. The region 
of 7 considered is not properly a restriction since we are 
interested in the small -> limit. All the arguments can 
be easily extended to the Lorentzian signature, that will 
be discussed later. For each face / (the 2-cells) of the 
2-complex a there is an associated integer spin jf . Faces 
arc oriented and bounded by a cycle of edges e (the 1- 
cells). Each edge bounding a face has a source vertex 
s(e) and a target vertex t(e), where source/target is rela- 
tive to the orientation of the face. To each edge associate 
SU (2) elements n e f (/ runs over the faces meeting at the 
edge e) and two source/target 

Spin{4) ~ SU(2) x 51/(2) ~ S*0(4) (5) 



fined as 




The sum is over spins j/ £ N, and the integrals are over 
the Spin{A) gauge variables and SU(2) variables labeling 
the edges. We dropped a face normalization factor which 
is usually taken as the dimension d(jf) = 2jf + 1 of the 
5J7(2)-irrcducible Hilbcrt space. Other normalizations 
are possible but irrelevant in the present analysis. 
The face amplitude Pj is given by 

Pf=^l[eefPef (7) 

where II denotes the ordered product of operators (ac- 
cording to the cycle of edges) and 

Pef = g e , s (e)Y[j fl n ef )(j f ,n ef \Y < <g-l ier (8) 

Here \j, n) is the SU(2) Bloch coherent state [25] for an- 
gular momentum 1 along the direction of n. The map Y 
gives the embedding of the SU{2) irreducible representa- 
tion space with spin j into the SU(2) x SU(2) represen- 
tation space with spins (i + ,i~), where 2 



In the canonical basis the embedding map Y is given by 
the Clcbsh-Gordan coefficient 

m + m~ \Y\j; m) = C(j + j~j:m + m~m). (10) 

Using the decomposition of the 50(4) variables g in self- 
dual and antiselfdual rotations (g + ,g~) € SU (2) x SU{2) 
and the factorization properties of SU(2) coherent states 
[26] we can write 

Pef = Ptf ® Pef (11) 

with 

P f = lU^, (12) 

and 

P?f = 9tie)Kf) mj hn e fr j h9%(e))- 1 (13) 



gauge group variables g e ,s(e)i 9e,t(e)- The variables n e f 
can also be interpreted as unit vectors n e / in 1R 3 , up to 
a phase ambiguity, by saying that n is a rotation that 
brings a reference direction to the direction of n. 

The spinfoam amplitude, or partition function, for the 
2-complcx a in the Bloch coherent state basis [24] is de- 



1 We introduced the notation [j, n) for the standard antilinear 
map e applied to \j,n). In the standard basis, the one that 
diagonalizcs J 2 and J z , it is given by the symbol ■ , e mm i = 

2 Because of equation (9) 7 is quantized to be rational. This re- 
striction is not present in the Lorentzian version of the theory. 
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where the Bloch coherent states are now in the fundamen- 
tal j = i representation. Define also the face amplitudes 
in the fundamental representation as 



/ 



P 



Jf-- 



f 



(14) 



a notation will be useful later. The partition function (6) 
can then be written in the form of a path integral 



W{a) 



E 

0/} 



dg ev / dn e f e 



with complex action given by 



5 = E 5 / = E lnP /- 



(15) 



(16) 



This is completely equivalent to (6). If we want to focus 
on the nondegenerate part of the state sum, we can elim- 
inate a small integration region in the neighborhood of a 
zero measure set defined as follows 



(defined on the infinite 2-complex) as the result of the 
iteration of some kind of renormalization group transfor- 
mation at the level of the spinfoam 'lattice'. The latter 
possibility is intriguing but, to our knowledge, a little far 
away from the present computational techniques. 

In order to see the effect of (17) on the partition func- 
tion, let us restrict our attention to a 2-complcx which 
is dual to a simplicial triangulation. Now every vertex is 
surrounded by five edges and ten faces, and every edge 
bounds four faces. Vertices are dual to 4-simplices, edges 
arc dual to tetrahedra, and faces are dual to triangles. 
The analysis parallels the one of [23, 31] at fixed, large 
spins, and will exclude the degenerate part of the state 
sum, as done in [23]. 

We are interested in making explicit the dependence on 
the Immirzi parameter, so let us decompose the action 
(16) using (9) in the following way 

S f = T j f (In P+ - In P- ) + j j (In P+ + In Pf) (18) 

or, equivalently, 



THE CONTINUUM 7 ->• LIMIT 

The continuum limit of the theory is defined as the 
infinite refinement of the 2-complex [11] (possibly un- 
dergoing a second order phase transition to a smooth 
spacetime), where we expect a large portion of spacetime 
to be described as the union of many small elementary 
regions (quanta) labeled by spins of order one. Differ- 
ently, here we define a continuum limit that is suitable 
to describe gravitational physics with a truncation of the 
theory on a finite cellular structure of spacetime, namely 
with an approximation of the full theory. The spinfoam 
amplitudes associated to finite graphs can be viewed as 
effective amplitudes obtained after a coarse-graining pro- 
cedure [27, 28] applied to the complete theory. But since 
the truncated amplitudes are not fundamental, there is 
no reason to keep fixed the Immirzi parameter to its 
'bare' value 70. The possibility of a renormalization of 
the Barbcro-Immirzi parameter has been recently advo- 
cated in different contexts [29, 30]. 

Here we explore the possibility of a running towards 
zero, 7 — > 0, simultaneously with the large spin regime 
of the theory. Thus we consider 



j — > 00, 7 — > 0, J7 = const. 



(17) 



and 7J is the macroscopic physical area in Planck area 
units. Notice that the Immirzi parameter controls the 
spacetime discreteness. In particular it controls the area 
gap and the spacing between area eigenvalues, thus the 
limit (17) is the continuum limit for the area operator. 
The correctness of this approximation has to be checked 
against concrete computations in specific examples, or 
possibly be justified and derived from the full amplitude 



S f = a f (lnP+ - lnPp + -a/(lnP+ + hxPf) (19) 

where we have defined the area of the triangle dual to / 

as 



«/ = 73 f, 



(20) 



according to the area spectrum of loop quantum gravity. 

We are interested in the evaluation of the amplitude 
(15) in a region of macroscopic areas af and in the limit 
7 — > with af fixed (so jf — > 00). Collecting all the face 
terms (19), let us write the full action as 



s = s° + -S' 

7 



and the path integral 



W{a) 



E 



dg ve / dn e f 



-s' 



(21) 



(22) 



where now the action S is viewed as a function of the ar- 
eas, and consequently the summation is over 7-multiplcs 
of the integer spins. The continuum approximation of 
this sum can be taken only if the sum is dominated by 
a region of areas that contains a sufficiently large num- 
ber of area eigenvalues. This condition can be satisfied 
only when we consider the more general, and more physi- 
cal, transition amplitudes associated to a 2-complex with 
boundary. This will be discussed later. For the moment, 
let us suppose we can do this approximation and write 



W(a) 



dcif j dg 



dn 



ef ' 



(23) 
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where we have dropped 3 one global factor 7 -1 per each 
face, that takes into account the measure of the small 
area intervals in the Ricmann sum. 

For the stationary phase evaluation of (23) we will have 
to take variations of the second term S' of the full action 
(21), because it is the one proportional to the large pa- 
rameter - — > 00. Being 5" complex, in addition to sta- 
tionary points we have also to find out the points that 
maximize the real part. The action S' has nonpositivc 
real part ReS" < 0, so the main contribution to the inte- 
gral over a/, g ev and n e f comes from the critical points 
where 

ReS" = 0. (24) 

Since both and PJ have nonpositive real part, the 
condition (24) holds for ReP^ = RePf = 0. Using that 
the scalar product of coherent states satisfies |(ni|ri2)| = 
i+ni-n-2 ^ wc g e ^. eas jjy that the cr itical points are the 
solutions to the following system 

9tv™ef = -gt v n e >f, (25) 
9e V nef = -g~, v n elf , (26) 

where e, e' arc adjacent edges in the face /, sharing the 
vertex v. In the last equations g ± acts in the vector repre- 
sentation on R 3 . If there are no solutions, the amplitude 
is exponentially suppressed. Using (25, 26) we have also 
that the brackets 

(n ef \(gte)- 1 9tK'f) = e< (27) 

reduce to simple phases, on the critical points (on-shcll). 

Furthermore, we must require that on the critical 
points the action is stationary. Varying S' with respect to 
50(4) group variables, and evaluating at a critical point 
we obtain the condition 

^S'\ crit =0 — ► J2*ftief = (28) 

fee 

which expresses the closure relation for the tetrahedron 
dual to the edge e. Variation with respect to the unit 
vectors n e f does not give further information, because it 
is automatically satisfied. Before completing the contin- 
uum approximation, let us open a short discussion on the 
geometrical interpretation of the critical points. 

The existence of critical points is related to the exis- 
tence of a triangulation where the areas 4 of triangles are 



3 We are only interested in the oscillatory behaviour of the integral. 

4 It is well known that a generic assignment of areas does not 
correspond to any Regge triangulation [32-34], In other words, 
in general there exists no assignment of lengths Zj, for the sides b 
of triangles such that aj = aj{l^,). The critical equations select 
only the Regge-like area configurations. 



specified by the set a/, and the unit normals to the trian- 
gles (viewed in the 3-dimcnsional space of each tetrahe- 
dron e), are specified by the set of unit vectors n e f. Thus 
the critical points have the nice interpretation of a 4- 
dimensional Regge manifold, that is a manifold endowed 
with a continuous, piecewise flat metric where curvature 
is distributional and concentrated on triangles. This can 
be seen as follows. First, define the unit vectors 

n vf = dtv^ef (29) 

where g^ v and n e f are a solution of the critical system 
(26) and the closure relation. Then the vectors 

Jff = jfffif (30) 

can be interpreted as the selfdual (+) and antisclfdual 
(— ) components of the Lie algebra valued discrete field 

Jvf = B v f + -B*f e so(4) (31) 

from which we can extract the field B v f. The field B v f 
codes the spacetime metric degrees of freedom because 

A v f = e 3 B* f = E vfbl A E vfb2 (32) 

is the simple area bivector 5 of the triangle /, in the frame 
of the 4-simplex v, where we have restored dimensional 
units. It is the simplicial version of e A e, and it is con- 
structed via the vectors £ R 4 along two sides, la- 
beled with b, of the triangle dual to /, in a flat coordinate 
patch covering the 4-simplex v. 

So far we have reconstructed the geometry of the 4- 
simplices out of a given solution (a/ , g ev ,Ti e f) of the crit- 
ical system (25, 26) and the closure relation (28) (see [31] 
for the 4-simplex reconstruction theorem) . The next step 
is to check if 4-simplices are consistently glued together, 
by identification of the shared tetrahedra, so that they 
form a proper spacetime triangulation, and (32) provides 
a consistent discrete tetrad over the simplicial manifold. 
By construction, the Euclidean geometry of the tetrahe- 
dron dual to e is the same when viewed from the source 
4-simplex s(e) or from the target 4-simplex t{e). 

However there is a discrete ambiguity in the way 4- 
simplices are glued [35]. Indeed from any given solution 
(a,f,g ev ,n e f) we can generate the finite set of parity- 
related solutions by a flip in the selfdual and antisclfdual 
group variables , g~ v in an arbitrary subset of vertices 
v. This is equivalent to a parity transformation, or reori- 
entation, of a subset of 4-simplices. Thus {af,g ev ,n e f) 



5 We speak about bivectors or elements in the Lie algebra as if they 
were the same thing, using the algebra isomorphism A 2 R 4 ~ 
Lie(0(4)). 
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labels a class of solutions, and the elements of this class 
arc specified with the set e v of orientations of 4-simplices. 
Now we have that two 4-simplices are glued consistently 
with a Regge triangulation only if their orientations in- 
duce opposite orientations in the common tetrahedron. If 
the last condition is true for every couple of vertices, the 
solution has the correct global orientation and describes 
a Regge manifold. In all the other cases it describes only 
a generalized Regge manifold, with no simple geometri- 
cal interpretation, and will produce a generalized Regge 
action, as we shall see. 

The spin connection is coded in the SO(A) holonomy 
U e that parallel transports the bivectors from one 4- 
simplcx to another along the edge e. It is the unique 
solution of the simplicial Cartan structure equation 



E, 



t(e),fb — U e E s ( e )j b 



(33) 



and defined 6 in terms of the on-shell gauge variables in 
a simple way: 



U e - 9e,t(e)9 eta ( e y 



Observe that (27) implies that the full holonomy around 
the loop (boundary of the face) has the form 



(36) 



where the first holonomy in the product is along the edge 
e such that s(e) = v, so the holonomy is based at v. Thus 
Uf is a SO(4) rotation that decomposes into a rotation 
along the plane of B v f, and a rotation along its dual 
orthogonal plane, the one of £?*y, and preserves both 
B v f and B*^ (bivectors have no definite location in their 
plane). The angles are given respectively by 



vef 



(37) 
(38) 



However, the Cartan structure equation (33) puts a fur- 
ther constraint on Uf. Apart from preserving B*^ it 
must preserve vectors in the plane of B*j. It follows that 



From its definition, it is evident that U e parallel transports the 
area bivectors 



U e A, 



(<=),/■ 



(34) 



Considering a single face, the previous equation still does not 
determine whether it parallel transports the side vectors E s ( e j j b 
as well, because we can do an arbitrary rotation of the target 
Et(e),fb along the plane of A t ^ j without affecting the condition 
(34). However (34) holds for any of the four faces bounded by 
e. It follows that U e parallel transports the side vectors as well, 
thus U e is the spin connection. 



<d*f vanishes. Interestingly, 9y is an angle with the inter- 
pretation of a torsion degree of freedom and accordingly 
it must vanish since the on-shell discrete connection U e 
is torsion- free (see also [23] for more details on torsion). 

It is finally time to get back to the small 7 expan- 
sion. The last thing to do is to take variations of S' with 
respect to the areas a/, but this does not give further 
restrictions, because 



dS' 
dcif 
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(39) 



is automatically satisfied on the critical points. Thus in 
the 7 — > limit the total action S must be evaluated at 
a critical point of 5", and using (27) we have 



S'|, 



S f 



f 



J- 



(40) 



(35) 6 



The angle 9/ is the generalized deficit angle, given by (2n 
minus) the sum of all the 4-dimensional dihedral angles 
'^f — 9~j between the neighboring tetrahedra that share 
the triangle dual to the face /. Notice that if the solution 
is such that the 4-simplices around a face are glued with 
consistent orientations then 6/ is the usual deficit angle 
of Regge calculus. 

The deficit angle codes the spacetime curvature con- 
centrated on the triangle. Its relation with the holonomy 
around the face is given by the previous formula (36) 
with 9} = 0. 

The quantity (40) is the Regge form [36, 37] of the 
action Sr for general relativity. However here what we 
get is a generalized Regge action, in the sense of Barrett 
and Foxon [38] , a situation similar to the Ponzano- Regge 
model for spinfoam gravity in three spacetime dimensions 
[39]. The exponential of (40) is the contribution to the 
partition function associated to a given critical point. In 
order to write the asymptotic approximation of the in- 
tegral (23) we have to sum over all critical points. In 
general this will take the form of a continuous sum over 
the set C of critical points (critical manifold) , with some 
measure /i that in principle can be computed by standard 
asymptotic analysis tools, in the following way. 

Suppose we are given an n-dimensional integral of the 
form 



I = Jd X a( X )e 



A/Or) 



(41) 



and we want to compute its asymptotic expansion for 
large positive A. Suppose the action / has nonpositive 
real part and a continuous set C (hypersurface) of critical 
points y G C where Ref(y) = and df/dx l \ y = 0. Then 
we have the following expansion: 



I 



C y ^ y '^m 



(y) 



(1 + 0(1/A)) (42) 
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where H is the Hessian matrix of / restricted to the 
directions normal to the critical surface with respect to 
some metric, and fi is the measure induced on the critical 
surface by the same metric. In the factor on the left of 
the integral sign, the action is evaluated on an arbitrary 
reference critical point tjq (but it is of course independent 
of this choice). 

Using the general formula (42) for / = S',a = exp(S' ), 
and A = I/7 we find for the partition function 

W(a) = I d^a f ,n ef ,g ev )e^ +0(j) (43) 
Jc 

where 0(7) denotes the 7-corrections to the amplitude 
that come from the next-to-leading orders of the asymp- 
totic approximation. ' The integration region C contains 
also the contributions from all the possible orientations 
e v of 4-simplices. 

THE SEMICLASSICAL H -t LIMIT 

More interestingly, let us parametrize the previous in- 
tegration (43) using length variables. Given that a class 
of parity-related critical points (af,n e f,g ev ) corresponds 
to a Regge triangulation 8 , there exists an assignment of 
lengths lb to the sides b of the triangles such that the 
areas a/ coincide with the areas computed out of the 
lengths h- 

aj=a f {l b ). (44) 

Exploiting this fact, we can parametrize the critical man- 
ifold C with the set of side lengths lb and with the joint 
orientations e = (ej., £2, ■ • ■) of the 4-simplices. Restoring 
the h dependence, we rewrite (43) as 

W(a) =J2[ d ^ e* S ' R(h) + 0(7), (45) 

and the generalized Regge action is now an explicit func- 
tion of the lengths 

ShW^afiWfQt) (46) 
/ 

as in the original classical formulation [36]. We have 
also made the sum over the orientations e explicit, so the 
generalized deficit angle at a face depends also on the 
particular orientation of the 4-simplices around the face. 



The corrections in 7 should contain also the error we make by 
approximating the sums with integrals. This is expected to decay 
at least exponentially. 

Only the consistently oriented elements of this class have to be 
interpreted as a proper Regge triangulation. 



The last expressions are a good starting point for tak- 
ing the semiclassical limit. Before discussing this, let us 
make a few comments on the effect of sending the Im- 
mirzi parameter to zero. The remarkable consequence of 
the continuum limit 7 — > is that the spinfoam ampli- 
tude reduces effectively, that is up to 7-corrections, to a 
quantization of Regge gravity [40] given by formula (45) 
where the fundamental variables are continuous lengths, 
the major difference being the nonzero contribution of the 
non-geometric orientations. The result resonates with 
the recent findings in the computation of the graviton 
propagator within loop quantum gravity [16-18]. In par- 
ticular, as shown in [18], the leading order (in the h ex- 
pansion) graviton propagator G(x, y) presents the same 
kind of 7-corrections, 

R + iX + -i 2 Y 
G(x, y) = n y V + h-corr. (47) 
\x - y\ l 

and only in the limit 7 — > the tensorial structure of the 
2-point function matches with the matrix of correlations 
R computed in quantum Regge gravity, and, even more 
interestingly, with the one given by standard pcrturba- 
tive gravity on flat space. In retrospect, at the light of 
the present general analysis the previous result (47) is 
much more clear. 

Suppose now we are interested in the semiclassical ex- 
pansion of the partition function (45). This corresponds 
to looking at a region of the integration domain where 
the areas are macroscopic, that is large as compared to 
the Planck area, 

a 4r » 1, (48) 
Ip 

or equivalcntly to the standard formal WKB expansion 
h — > 0. Remember that we are still in the small Immirzi 
parameter regime. We can state the regime in a more sug- 
gestive fashion by introducing two physical scales. One 
is the length scale of quantum gravity, identified with the 
Planck length 

Iqg = lp, (49) 

the other is the scale of loop quantum geometry, that is 
the scale where we can 'see' the discreteness of spacetime 

l QG = yfflp. (50) 

Thus the regime of the spinfoam path-integral we are 
looking for is expressed by the following relation: 

I > he > Ilqg (51) 

where I is the typical linear scale of each 4-simplcx. 
This regime can be selected by appropriate semiclassi- 
cal boundary conditions in the transition amplitudes, as 
explained later. 
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The classical equations of motion are obtained by vary- 
ing (46) with respect to the lengths. Using also the 
Schlafli identity [36] , which tells that the variation of the 
deficit angles does not contribute to the total variation 
of the action, these are the well known Regge equations 

£^=0, (52) 

a discrete version of the continuum Einstein equations in 
vacuum 

JV = °' ( 53 ) 

that is of the vanishing of the Ricci tensor. Notice that 
the interpretation of (52) as Einstein equations is clearer 
(probably possible only) when e is a consistent global 
orientation. 

The Regge equations (52) enforce a relation between 
the generalized deficit angles of different faces. The in- 
tegral (45) is dominated by its stationary 'trajectories', 
namely by the sets of lengths lb which are a solution of 
the Regge equations (52). However, in order to pick up 
a single classical trajectory we need to specify appro- 
priately the boundary conditions, or in other words we 
have to consider the transition amplitudes. The spin- 
foam boundary formalism for the transition amplitudes 
is briefly reviewed after the following Lorentzian section. 

LORENTZIAN SIGNATURE 

The model in its Lorentian version [7] has pretty much 
the same form as the Euclidean one, mutatis mutandis. 
The gauge group is now 

SX(2,C) ~ 50(1,3). (54) 

As before, the partition function is 




with face amplitude Pf given by the trace of operators 

Pf^^eefPef (56) 

and where the elementary edge operator is the infinite 
dimensional matrix 

Pel = g e ,s(e)Y[jf,nef){jf,n e f\Ylgl itie) . (57) 

Let us explain the previous formulas. As in the Eu- 
clidean theory, we have the same SU(2) spins jf an d 
Bloch SU{2) coherent states labeled by n e f. The gauge 
group variables g ev are elements in SL{2 1 (D), integrated 
with the Haar measure, and the redundant integrations 



(infinite volumes of SL(2, (D)) have to be dropped in or- 
der to make the amplitude potentially finite [41]. The 
map Y is constructed through the principal series of uni- 
tary representations Ti^'^ of SL(2, (D), which are infinite 
dimensional and labeled by an integer k and a positive 
real number p. We have the useful orthogonal decompo- 
sition 

W (w, =0^ (58) 
k>j 

where the subspace T-Lj carries a representation of spin 
j of the SU{2) subgroup of SL{2, (D) rotations that leave 
a fixed timelikc 4-vector invariant. The map Y in (57) is 
defined in terms of this decomposition as the isometric 
embedding of SU (2) irreducible representations j in the 
lowest weight term in (58), namely k = j. Moreover, the 
Lorentzian model imposes the following constraint on the 
labels k and p: 

P = lk, (59) 

which is the analogous of (9). Notice that the previ- 
ous equation does not impose a quantization of the Im- 
mirzi parameter, that can be any positive number, a 
comforting feature from the perspective of the contact 
with Hamiltonian loop quantum gravity, where 7 is not 
quantized as well. Hence the map Y is the embedding 
of the 5(7(2) representations j into the SL(2, (D) repre- 
sentations (j, 7j). This completes the definition of the 
Lorentzian model. 

In order to write a path integral form of the model, it 
is useful to realize the Hilbert space T-[( k ' p ) in the usual 
way as the space of homogeneous functions of z G (D 2 
with a specific measure f2 that defines the scalar product 

(/,«?) = f dzn{z)J(z)g(z) (60) 

and the integration on the complex projective line (DP 1 
makes sense because f2 is such that the integrand 1-form 
is invariant under rescaling of z. Thus z = (z\,Z2) € 
(D 2 can be interpreted as homogeneous coordinates [z 1 : 
Z2] for (DP 1 . It is a simple exercise in representation 
theory to compute explicitely the integral kernel of the 
edge operator (57), which has the following expression 

where the function Q e f is slowly varying. Now we define 
the analogous of (19) as 

5 e/ =ia / logP e ° / + ^logP e ' / (62) 
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with 



P. 



lfL(e) Z s(e),/l 



e f | t I 

Z s {e)j\9e,s(e) K/) (n e f W et[e) \z t (e),f) 



\9l. s(e fs(e)j\\9l Ae ft(e)j\ 

and the total action for the 2-complcx a as 



EE*-/ 



S° + -S' 



(63) 
(64) 

(65) 



where in the last equality we have collected all the terms 
proportional to •=■. From the form of the edge kernel (61) 
it is clear that this action formulation is equivalent to the 
original definition of the partition function (55), and our 
final form is 

W(a) = j dg ev J dn e j J dz v juje s (66) 



with measure density given by 



(67) 



Observe that we have one integration dz v f on (DP 1 per 
each couple (v, /) of the 2-complcx. The action (65) and 
the variables used are a simple generalization to an ar- 
bitrary 2-complcx a of the single vertex action of [42]. 
Notice also that 5° is purely imaginary while S' is com- 
plex. While the previous definitions are general, from 
now on we restrict our attention to the particular case of 
a simplicial 2-complex. 

In the continuum limit we look for the critical points 
of the complex action 5". Imposing ReS" = and the 
vanishing of variations of S' with respect to z v f we find 
two equations 



g e v [n e f) 



\glv z vf\ 
\glv z vf\ 

\gt'v z vf\ j-i 

\gevZyf | 



9e'v\n e >j) x phase 



-9e 



n e ' f) x phase 



(68) 
(69) 



for the critical points, where e and e' are two adjacent 
edges in the boundary of the face /, both bounded by 
the vertex v. The requirement of stationarity of S' with 
respect to the n e f variables gives the trivial identity n e / — 
n e f = 0. The variation of S' with respect to the areas 
a/ is automatically zero on the critical points, because 
of the torsionless condition, similarly to the Euclidean 
case. Varying with respect to g ev and using the previous 
critical equations we get the closure condition in each 
tetrahedron: 



E a /™ e / 

fee 



0. 



(70) 



The solutions to those equations can be interpreted as 
4-dimcntional Lorcntzian Rcggc manifolds, using the fol- 
lowing construction. Given a solution, consider g ev , n e f 
and out of them define two null (lightlike) 4-vectors in 
Minkowski space 



N v f = g e v(l, fief), 

N v f = g eu (l, -fi e f). 



(71) 
(72) 



Taking their wedge product, define the spacclikc simple 
area bivectors 9 



A vf = cif(N v f A Nvf)* G so(l,3) 



(73) 



that code the geometry of a curved spacetime triangu- 
lated with flat 4-simplices. Indeed for each vertex v, 
these bivectors describe completely the geometry of a 
Lorentzian flat 4-simplex with spacelike boundary: there 
exists a Lorentzian 4-simplex with area bivectors defined 
as (73), whereas in the case the variables are not a solu- 
tion of the critical equations such a 4-simplex does not 
exist (see [42] , [43] on the reconstruction theorems for the 
4-simplcx). Two neighboring 4-simpliccs are described 
by two sets of ten bivectors A s r e \j and A t ^j. The 4- 
simplices are glued together and share the tetrahedron 
dual to e. The 3-dimensional Euclidean geometry of the 
common tetrahedron is the same when viewed from s(e) 
or t(e) by construction of the model; this is due to the 
fact that in the edge amplitude (57) the source and tar- 
get unit vectors are not independent but are both equal 
to n e f. However, as in the Euclidean signature we can 
glue two 4-simplices in two ways, one with the correct 
orientation, the other with the 'wrong' orientation, and 
both solutions are allowed. 

Thus (73) defines the metric of a (generalized) 
Lorentzian Regge manifold, with Levi-Civita connection 
coded in the parallel transport 



U e 



5e,t(e)fl , e , s (e)' 



(74) 



interpreted as the path-ordered exponential of the 
Lorentz connection from one 4-simplex to another. The 
boundary of every 4-simplex is formed by five spacelike 
tetrahedra. It can be shown [42] that the real ratio 



CXp Uyf 



(75) 



of equations (68), (69) yields the Lorcntzian dihedral 'an- 
gle' 6 v f between the two tetrahedra dual to e and e' in 
the boundary of the 4-simplex v. It is not an angle, but 
rather a boost parameter. Its hyperbolic cosine is the 



9 We always trade bivectors for Lie algebra elements and viceversa, 
using the isomorphism A 2 R 1,3 ~ Lie(0(l, 3)). 
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scalar product of the unit timelike outward normals to 
the tetrahcdra, computed in the 4-simplex v, up to a 
sign. If one normal is future pointing and the other past 
pointing the dihedral angle is defined to be positive. If 
both are future pointing, or past pointing, the angle is 
defined negative. 

The second term S' vanishes on the critical points, 
because of the torsionlcss condition (analogous to the 
Euclidean case). So the total action gives the on-shell 
contribution 

S\cnt. = S R = i} j a/0/ (76) 
/ 

where we have defined 

which is the generalized Lorentzian deficit angle at the 
face /, dual to a spacelike triangle. Notice that similarly 
to the Euclidean signature, also in a Lorentzian triangu- 
lation the deficit angle is zero for a flat connection. The 
on-shell action (76) is the (generalized) Regge action for 
Lorentzian gravity. It is worth observing again that the 
critical equations admit solutions only for Regge- like con- 
figurations of the areas, because we cannot reconstruct 
the geometry (73) using areas which are not admissible. 
From here, the analysis of the continuum and semiclas- 
sical limits proceeds exactly as in the case of Euclidean 
signature. 



GENERAL BOUNDARY AMPLITUDES 

The predictions of the spinfoam theory are the transi- 
tion amplitudes. These are associated to general bound- 
aries of the 2-complexes [44]. Given a 2-complex with 
boundary (the boundary is interpreted to be spacclike in 
the Lorentzian theory), its boundary graph T is an ab- 
stract oriented graph made of links I (where the external 
faces end) and nodes n (where the external edges end). 
The boundary graph inherits its labeling and orienta- 
tion from the external faces and edges. The set formed 
by the boundary graph T, the spins ji associated to the 
links, and unit vectors n n i associated to the nodes is the 
boundary data. Notice that for a simplicial 2-complcx 
the boundary graph contains only 4-valent nodes. So for 
each node we have exactly four labels n n \ , where I labels 
the four links bounded by the node n. These are inter- 
preted as the four normals to the faces of the tetrahedron 
dual to the node. 

The spinfoam transition amplitude for the 2-complex a 
with boundary T in the Bloch coherent state basis [24] is 
a functional of the boundary data given by the following 



expression 




where now the face amplitude is split in a boundary (rel- 
ative to external faces I), and a bulk part (relative to 
internal faces /), and accordingly the two products are 
over the external and internal faces respectively. Let us 
explain better the last formula. First, the external face 
amplitude Pi is the same as (7) except that for the edges 
ending on the boundary we have only 'half of (8), namely 

P ef = (j f ,n ef \Y^gl Me) (79) 

or 

Pef = 9 e ,s( e )Y\jf,n e f) (80) 

depending on the orientation of the edge induced by the 
face orientation. Second, the summation is over the sole 
internal spins jf € N, and the integrals are over the 
Spin(A) gauge variables in the Euclidean, or SL(2, (D) in 
the Lorentzian, and also over SU (2) variables (the unit 
vectors) labeling the internal edges. 

The amplitude (78) for a 2-complex with boundary has 
the nice interpretation as transition amplitude associated 
to a spin- network supported on the boundary graph T. 
The canonical loop quantum gravity Hilbert space asso- 
ciated to the graph T is 

H r = L 2 (SU(2) L /SU(2) N ) (81) 

with L the number of links and N the number of nodes. 
If r has two components T = Ti n U T ouf we could in- 
terpret the transition amplitude (78) as a standard tran- 
sition from an initial state in Hr in to a final state in 
%T OU f However it could be more correct to call those 
amplitudes extraction amplitudes, because in the spin- 
foam paradigm they are supposed to 'extract' the physi- 
cal states by projecting the kincmatical boundary states 
down to the Hilbert space of the solutions of the Hamil- 
tonian constraint. 

In the state space (81) there exists an overcomplete 
basis of spin- networks, the Bloch coherent state basis in- 
troduced in quantum gravity by Livinc and Spcziale [24] 
labeled precisely by ji and n n i- Thus (78) can be in- 
terpreted as the covariant transition amplitude for the 
quantum states of 3-geometry given by the loop quan- 
tum gravity states \ji, n n i), 

W(j h n ne ) = {W\ji,n n i}. (82) 

The semiclassical analysis of the previous section is 
straightforwardly generalized to a 2-complex with bound- 
ary, namely to the transition amplitude (78), once we 
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write it in the path integral form. For example, in the 
Lorcntzian case we rewrite 

W(j h n ne ) = ^2 J d 9veJ dn ef J dz vf u>e^ Sl+s f. (83) 
{it } 

The boundary formalism allows to select the contin- 
uum/semiclassical regime in the following way. First, 
with a little abuse of notation, write the boundary func- 
tional in terms of the areas as 

W{a u n ne ). (84) 

Then we look at the behavior of this transition amplitude 
when all the boundary areas ai are macroscopic, in the 
limit 7 — > with ai fixed. This is the 

a > Iqg > 1 Iqg (85) 

regime. It is formally the simultaneous semiclassical limit 
h — > and continuum limit 7 — > 0, and corresponds to 
looking at macroscopic geometries and at the same time 
neglecting the quantum discreteness of spacetime 10 . Sup- 
pose also that the path integral is dominated by areas 
of the same order of magnitude of the boundary areas: 
ai ~ a/. In other words, the macroscopic boundary state 
must enforce macroscopic areas in the bulk of the trian- 
gulation, a condition that can be checked for the specific 
geometry chosen. 

If we go through the previous analysis of the partition 
function, it becomes clear that the boundary amplitude 
receives the dominant contribution from the solutions of 
the generalized Reggc equations (52) compatible with the 
specified boundary data (see [45, 46] on the initial value 
problem). The nature and the number of solutions has 
to be studied in detail case by case. As an example, in 
the case there is only one solution of the Reggc equations 
the behaviour of the boundary amplitude is 11 

W(ai,n ne ) ~ Ae* SRiai ' nne) + s e -£ s «( a <.™-) (86) 

where now 5^ is the Hamilton function, that is the ac- 
tion evaluated at the classical trajectory determined by 
the boundary data via Einstein equations. The two expo- 
nentials with reversed global sign are related to the two 
opposite global orientations of the triangulation, and the 
quantum theory sums over them with a different weight, 
even if classically they represent the same un-oriented 
spacetime. This result for general triangulations is sim- 
ilar to the single vertex results of [31, 42, 47, 48]. In 



Notice that also in the h — > limit with area fixed the spectrum 
becomes quasi-continuous, thus our separation of the semiclassi- 
cal limit and the continuous spectrum limit is artificial and only 
a matter of speaking. 

The prefactors A and B are non-oscillatory and slowly varying. 



general we expect non-zero contributions also from the 
non-geometric orientations that must be added to (86); 
this may happen if there are solutions of the generalized 
Regge equations with the same boundary data. 

Notice that the boundary data a / , n n i may fail to be a 
consistent labeling of the boundary spin- network, namely 
there may not exist a 3-dimcnsional (spacelike, in the 
Lorentzian case) triangulation with the prescribed ar- 
eas and unit vectors, the latter interpreted as the unit 
normals to the faces of the tetrahedra. Due to the ge- 
ometric interpretation of the critical points, in this case 
the transition amplitude is exponentially suppressed. On 
the other hand, a consistent set of boundary data can be 
equivalently mapped to the set of lengths of the boundary 
triangulation. This set of lengths is a Dirichelet bound- 
ary condition for the Regge equations of motion to be 
used to determine the classical solution in the interior 
and evaluate the Hamilton function in (86). 

The result (86) is what we expected from a theory of 
quantum gravity, and is the concrete realization of the 
equation (4) in the introduction. 

CONCLUSIONS AND OUTLOOK 

It this paper we have discussed a proposal for the 
semiclassical limit of Euclidean and Lorcntzian spin- 
foams truncated to an arbitrary, finite triangulation (2- 
complcx), where most calculations are done. We find 
(equation (86)) that the transition amplitudes yield the 
exponential of the Hamilton function of Rcggc-Einstein 
general relativity, as expected, up to /i-corrections and 7- 
corrections in the simultaneous semiclassical and contin- 
uum limit. The first corrections correspond to the stan- 
dard WKB expansion of the path integral. The latter are 
new and are the effect of the discreteness of geometry, in 
the sense that the spectra of areas and volumes are dis- 
crete and the discreteness is controlled by the Immirzi 
parameter 7. The continuum limit we take is a looser 
concept of the 'full' continuum limit defined as the phase 
transition to a smooth spacetime manifold. However, 
working with a fixed, finite triangulation the only way 
of taking a continuum limit is to look at the continuous 
spectrum limit of the fundamental geometric operators. 
We have done this by letting 7 run to zero, keeping fixed 
the macroscopic areas a/ (this is still a tentative pro- 
posal and must be further investigated). Remarkably, as 
explained in the paper, the result sheds new light on the 
previous calculation of the graviton propagator with the 
'new models' [19], where the same kind of 7-corrections 
to the standard perturbative tree-level propagator have 
been found (one could speculate on potentially observ- 
able signatures of those pure LQG corrections). 

It is also very interesting to remark that essentially 
the same continuum limit was considered by Bojowald 
[49] in the context of loop quantum cosmology. Quoting 
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its abstract: "standard quantum cosmology is shown to 
be the sim.ultanpnv.s lim.it 7 — >■ j — »■ no »/ /on/) quantum. 



[s 



COSlhnlogtfl , a strong analogy it is worth studying further. 
Wc did not study in great detail the physical effects of 



[9 



i) tl |LC contributions with inconsistent orientations, and 

ii) the degenerate triangulations. In principle, these do 
not pose a real problem: for example, see [50] on the 
fate of the sign-reversed exponentials in the bulk of the 
triangulation, and how a coherent boundary state peaked 
on the appropriate intrinsic and extrinsic curvature could 
be able to select the correct global orientation (see also 
[51]). In the Lorentzian setting, the issue of the wrong 
orientations has been intimately related to the issue of 
implementing causality in the spinfoam models [52-54]. 

Finally, we have also disregarded the potential 'in- 
frared' divergencies associated to bubbles in the foam, 
similar to the loops of pcrturbative QFT, for which a 
suitable regularization and renormalization scheme [55- 
58] is required. 
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